Electron transport in parallel quantum wires with random potentials 
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We study an electron transport property in two parallel quantum wires with random potentials. 
Assuming the same microscopic parameters for both wires, we focus on the relationship between 
inter-wire interaction and electron backward scattering by random potentials at low energy regime. 
Our analytical and numerical calculations show that the Drude weight, a measure of the electron 
transport, is influenced by inter- wire interaction and random potential independently, and little 
coupling between those two is observed, which is in contrast to a deep relationship between up- and 
down-spin interactions and random potentials in a single wire. It leads to that inter-wire interactions 
do not have a great influence on the Anderson localization in each wire. 



The problem of electronic transport has been one of the 
major issues in the history of condensed matter physics. 
There are always two key factors: electron-electron in- 
teraction and random potential. The importance of their 
roles on the electronic transport is extremely high in low 
dimensions because of large quantum fluctuations and 
strong localizations. 

In this paper we pick up a one-dimensional system of 
parallel quantum wires. If no hopping between the wires 
is allowed, electrons on one wire recognize the presence 
of the neighboring wire only through inter- wire coulomb 
interaction. 

In one-dimensional system, interacting electrons form 
so-called Luttinger liquid and the important theoretical 
works have been done on its transport properties. Kane 
and Fisher and Furusaki and Nagaosa studied Lut- 
tinger liquid placed in a single (and double) barrier sys- 
tem and showed the various interesting aspects of con- 
ductance. 

Luttinger liquid with random potential was studied, 
for example in Ref . [3-5] . Electrons interacting each other 
with repulsive force are shown all localized even by an in- 
finitesimally weak random potential 0. There is in fact 
an interplay between the interaction and the random po- 
tential, but the transport property would not reflect it 
because the system is always an insulator. The interplay 
becomes significant in considering the transport in finite 
size systems that could be metal Q ■ Reference [j| claims 
that repulsive interaction would suppress the effects of 
random potentials and hence enhance the electronic cur- 
rent. 

If we put two equivalent wires in a parallel position, 
another interesting phenomena is expected. When no 
electron transferring between the wires is assumed, the 
electrons on each wire are confined in the line, but the 
inter- wire interactions couple the electrons of both wires. 
Electronic current in one wire would drag the electrons in 



the other wire, which is called coulomb drag effect. For 
the parallel wires, some theoretical works have been done 
concerning the temperature and electron reservoir depen- 
dence of transresistance, transconductance and transsus- 
ceptibility for clean wires as well as wires with periodic 
potentials [6-8]. 

In this paper we focus our attentions on, if any, the 
interplay between inter-wire interaction and random po- 
tential. We would like to see how the electronic transport 
in a dirty wire is modified by the neighboring wire that is 
also dirty. For this purpose we calculated Drude weight 
of a wire in the presence of neighboring wire and try to 
seen how it varies by changing both the interaction and 
random potential, expecting it to provide us informa- 
tions on whether the inter-wire interaction has any effect 
on the Anderson localization within the wire. 

The model we consider here describes spinless electrons 
in two identical parallel one-dimensional wires with ran- 
dom potentials. It could be interpreted that we are treat- 
ing parallel wires in a strong magnetic field that polarizes 
all the spins. Inter-wire and intra-wire interactions be- 
tween electrons are also assumed. The Hamiltonian is 
then given by 



H = -t^2(cl tW Ci + i >w + h.c.) + U n h+ i 

i,w i 
i,w i,w 



(i) 



where t is the hopping integral, U is the inter-wire in- 
teraction, V is the intra-wire interaction, R w is random 
potential on the site i and w = ± presents the wire index. 

In the low energy regime, the kinetic-energy part be- 
comes 



Ho = itvfL 



r,w,k 



(2) 
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where r = ± is an index for the right/left moving elec- 
trons, v F = 2i sin(fcp) is the Fermi velocity, p TtW {k) is the 
electron density in the momentum space. The small mo- 
mentum transfer part of intra- and inter- wire interactions 
is given by 



Hf = L 1 i.925 w ,w' -w>)p+,w(k)p-, w >(-k) 

k,w,w' 

+ (2L)" 1 (946 w,w' w,—w' 



r.k.w.w' 



(3) 



where g 2 = 2V"(1 — cos(2fcp)) and 54 = 2V are the cou- 
pling constants of intra- wire interactions, and 92 = U 
and 54 = U/2 are of inter-wire interaction. The 2k F 
backward scattering between the wires is given by 



{x)ip_ tW (x), 



(4) 

where gi = U . The total Hamiltonian of clean wires is 
then given by H c = Ho + Hf + Tib, which is obviously 
equivalent to an isospin-1/2 single-channel system. 

Following the standard procedure of bozoniza- 
tion, we represent fermionic fields "4> r ,w(x) as 
i>r,w(x) = (2iray 1/2 exp[ir{k F x - <j> w {x)} + 
iO w (x)] by introducing bosonic fields 4> w (x) = 



»l9l/2 



[p+,w(q) + P-,w(q)] and their con- 



jugates IL w (x) = L- 1 J2 q e-^- a ^ 2 [p +tW {q) - p-, w (q)}, 
where a = +0 and d x 9 w {x) — i{n. w (x). In this 
bozon representation and the transformations 

</> c± = (0+ ± 0_)/V2,n c± = (n+ ± n_)/V2, the 

Hamiltonian becomes 



+ 



+ 
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J dx[K c _n 2 n 2 c _ - 
— T / dxcos(V8(f> c 
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Kr 
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where 



u c ± = v F [(l 



94, ±94 ^2 _ ^ 92 ± 92 j 2 jl/2 



K c± = 



2irv F 2nv F 
2irv F + ff 4 ± 94 - 92 T 92 ji/2 
2irv F + gi ± gi + g 2 ± 92 



(5) 

(6) 
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When the interactions are weak, u c ± and K c ± are given 

by 



U c ± 



K c ± 



v F (l 



V 



± 



u 



wv F 4ttv f 



V U 
1 (1 - cos(2fc F )) T t ■ 
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(8) 
(9) 



We now introduce quenched random potentials, which 
are assumed to be parameterized by two Gaussian ran- 
dom fields r\ w and £ w , where 77 is real and £ is complex. 



These two fields represent respectively the forward and 
backward electron scattterings from the random poten- 
tials, and have the Gaussian distributions: 

P Vw = cxp[-(2A / )- 1 J dx^Kx)} (10) 

P iw = cxp[-(2A b )- 1 J dx\U^)\ 2 }- (11) 

The scattering terms with the random potentials arc 
given by 

Hfrand = ^ / dx Vw(x)[p+, w (x) + P- iW (x)], (12) 

Hbrand = X] / ^ ^+,™ >™ ( X ) 

+ &(x)tf- >w (xW+, w (x)]. (13) 

It should be noted that introducing random potentials on 
each wire makes the present model variant from its resem- 
blance, the isospin-1/2 single-channel model. This is be- 
cause in the former case the profiles of the potentials that 
the electrons feel on each wire are different even though 
their statistical characteristics (average depth and stan- 
dard deviation) are the same. In the latter model, on 
the other hand, the electrons of up and down spins feel 
exactly the same potentials. Our two-wire model with 
random potentials actually corresponds to an isospin-1/2 
single-channel model with random potentials and a ran- 
dom magnetic field. 

Using the conventional replica trick, we integrate out 
the random field and obtain the following replica actions: 



S 
So 



So + Si + S2, 

n 1 

El drdx — 
. J 2tt 



(14) 



x (-!-[— (dr4>U(x,r)) 2 +u c+ (d x 4>l + (x,T)) 2 ] 

+ -^[—(d T fc_(^T)) 2 +U c _(dJl_(x,T)) 2 }), 

J\ C _ 1l C — 

(15) 

n „ _ 2 

ST= -]T drdr'dxdx'^-, 

i,3 

XCO S (V8[4>i_(x,T)-4>i_(x , ,T')]) 
-4A f (^f\x-x'\ 



x e 



(16) 



S 2 = -Y^ / drdr dx^-^ 

*,3 



x cos(V2[4>i_{x,r) - 4r> c _{x,T')}) 
x cos(V2[4>i + (x,T) - 4> 3 c+ {x,t')]). 



(17) 



where i and j are the replica indices. To derive the replica 
actions, we have shifted the bozon field as 



4> c + = 4> c + + vtot(x), 



(18) 
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Vtot(x) 



K, 



c+ 



V2u c 

K r - 



dzr] to t(z), 
dz7] Te i{z), 



y/2u c - 

Vtot(x) = r)+(x) + r]-{x), 
Vrd(x) = T]+(x) - -q-ix). 

For this action, we derived the following rcnormalization 
equations: 



dK c+ {l) 
dl 

du c+ (l) 



K, 



U 



-(K c - + K c+ )A b (l), 



c+ 



= K C+ (K C 

dl u c+ 



dKc-Q) 
dl 



K c+ )A b (l), 



2a 



(19) 
(20) 



-[(K c _+K c+ )A b (l) + —K c „gi(l)}, 
7T 



(21) 
(22) 



= —^——K C -(K C - + K c+ )A b (l) 7 
dl u c _ 

dA b {l) 



dl 
dgxjl) 
dl 



(3-K c _-K c+ )A b (l), 
2(1 (J), 



(23) 
(24) 



where /i = a/3n. 

We hereafter focus on the Drude weight D w of wire 
w(= ±), which describes the electronic transport within 
the wire +(— ) in the presence of the neighboring wire 
— (+). The Drude weight D± can be determined by 
the current-current correlation function (J±) = ((J c + ± 
J c _) 2 /2), where J c ± is symmetric/antisymmetric current 
densities. Since we assumed qualitative equivalence of 
both wires, we have D + = D- = D and therefore 



D = (( J 2 + ) + (J c 2 _))/2. 



(25) 



With no inter-wire interaction and random potential, we 
get from the quadratic form of the current d T 4> c ± in the 
action 



D = K c+ u cJ , 



K n -u n 



(26) 



omitting the irrelevant numerical factors. The qualitative 
change of D by the inclusion of the interaction and the 
random potential can be seen from the behaviors of K 
and u in the RG equations. 

At first, g\± always makes K c smaller as the renor- 
malization proceeds in Eq. I|21|) , indicating that the pres- 
ence of the inter- wire interaction gives the smaller Drude 
weight irrespective to the sign of the interaction. This is 
the dragging effect, which is stronger for larger interac- 
tion. 

Although the full dependence of Drude weight on the 
microscopic parameters has to be calculated by numeri- 
cally integrating the RG equations, we take a rough view 
on the behavior of the parameter A b neglecting the I de- 
pendence of the K on the r.h.s of Eq. H'-'il) . Integrating 



the equation from I = to I = ln(L/et) we get 



Ah 



A w (£/a) 1+2y[1 " cos(2fcF)1/ ™ F , 



(27) 



where we used Eq. (0 and A b Q is the initial value. Note 
that the inter-wire interaction U has no direct effect on 
the random potential parameter Af,, which means that 
the major role of U on the electronic transport is to drag 
the neighboring wire electrons. 

This is quite different from the case of the spin-full sin- 
gle wire model which is equivalent to the present model 
until the random potentials are introduced as stated 
above. The interaction dependence of the random po- 
tential parameter for the spin-full single wire is given by 
the following form (see Ref.||). 



A = A (L/a 



l-U / (ttv f )+2V[1-2 cos(2k F )]/iTVf 



(28) 



where U and V are the on-site and nearest-neighbor inter- 
actions, respectively. The important difference between 
Eq. (|2*7jl and Eq. l(2*5|l is the appearance of U in the 
latter. In the spin-full single wire, U, the interaction be- 
tween up and down spins, would effectively increase the 
random potential parameter A, which can be understood 
by the following argument: The spin- 1/2 electron system 
in one dimension has the SDW-dominated ground state 
and the SDW fluctuation becomes larger with increasing 
U, whereas it suppresses CDW fluctuations. Therefore 
the pinning effect due to impurities is weaker for larger 
U. Since the last term in the r.h.s. of Eq. I|21|) does not 
exist in the RG equations for the spin-full single wire case, 
meaning the absence of the dragging effect, the influence 
of U on the Drude weight D appears solely through A. 
Consequently U would suppress A and therefore D. In 
other words the repulsive interactions between up- and 
down- spins help the transport. The different roles of in- 
teractions on the localization in our model and the spin- 
full single-wire model can also be seen in the localization 
length £i oc . Following Ref. Q, we can estimate £i oc as 
follows, 



Cloc/a ' 



' (£j / a ) 1 - 2l/ [ 1 - cos ( 2fc F)]/™F 

(spinless two-wire model) 

^(°) / a y + U/(Trv F )-2V[l-cos(2k F )]/Tvv F 

(spinfull single wire model) 



(29) 



where a is the lattice constant and £^ is the localization 
length of non-interacting system. 

In order to check the above result numerically, we per- 
formed Monte Carlo calculations for the Hamiltonian QJ 
using the world-line algorithm. The random potential on 
each site is assumed to be a random number uniformly 
generated between —R and R. The simulation was per- 
formed in the zero winding number mode, based on the 
periodic boundary condition. While it is necessary to cal- 
culate the current-current correlation function to obtain 
Drude weight, its zero-frequency part < J(uj)J(—lu) > is 



4 



always zero because we worked in the zero winding num- 
ber space. We therefore extrapolated < J{m)J(—u>) > 
to uj — > limit and determined the Drude weight [lCj . 
The simulations are performed for wires of 20 sites and 8 
particles each. The random potential average were took 
over 100 samples to achieve satisfactory precision. 

Figure 1 shows the Drude weight as a function of the 
inter- wire interaction U, at the temperature (2/9)t. The 
intra- wire interaction V is fixed: V — (2/3)t. We see 
the monotonic reduction of D for either increasing U or 
increasing R. Let us normalize the Drude weight with 
the value of the clean system, i.e. D(U,R)/D(U,R = 0). 
Then we find the normalized D has little dependence 
on U (see the inset of Fig.l), which supports an idea 
that D(U, R) can be approximately written as a product 
of two functions, one of U and the other of R. If we 
write this as D(U, R) ~ F 1 (U)F 2 (R), F 1 and F 2 are both 
decreasing functions. The point here is that U and R are 
not coupled in the Drude weight, which is exactly what 
we have on the RG calculations. 

In conclusion we studied a transport property of spin- 
polarized electrons moving in parallel one-dimensional 



wires with random potentials. We in particular focused 
on the Drude weight, a measure of the transport in one 
wire under the influence of the other. Analytical and 
numerical calculations were done to see how the Drude 
weight varies as a function of the inter-wire interaction 
or of the random potential. We found both the inter- 
action and the potential reduce the Drude weight: the 
former is due to the coulomb drag and the latter to the 
localization. An interesting point is that there is no no- 
ticeable interplay between these two. This is in contrast 
to the case of full-spin electrons in a single wire, where 
interactions between up- and down-spin electrons would 
weaken the random potential effectively. Little effect of 
the inter-wire interaction on the strength of the random 
potential means that the property of Anderson localiza- 
tion in each wire has little dependence on the inter-wire 
interaction, although it is greatly affected by the intra- 
wire interaction and by the inter- wire hopping 9J. 
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FIG. 1: Drude weight D as a function of inter- wire interaction 
U for the spinless double wire system. The random potential 
varies between — R and R and t and V are fixed to be 1.5 and 
1 respectively. 



